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INTRODUCTION 


Cooperative watershed research between the Agricultural Research Service, 
U. S. Department of Agriculture, and the Bureau of Land Management, U. S. 
Department of the Interior, was initiated in 1968 under Cooperative 
Agreement No. 14-11-0001-4162(N). Also, the Memorandum of Understanding, 
dated July 6, 1961, which is a part of the Cooperative Agreement, specifies 
the overall responsibility of each agency. ‘ 


This interim report summarizes progress and results from January l 

through December 31, 1973, as outlined in the work plan for F.Y. 1974. 

The report also describes the proposed activities and changes in objectives 
for consideration in the F.Y. 1975 work plan. 


Significant changes in personnel at the Northwest Watershed Research 
Center were the addition of Dr. Clayton Hanson, Agricultural Engineer, 

to the professional staff, and the permanent appointment of several 
technicians to the support staff, formerly on temporary or part time 
appointments. These actions will provide greater competence in rangeland 
evapotranspiration studies and greater stability to the staff. 


The Reynolds Creek Watershed is included within an area in southwest 
Idaho, which has been classified as needing management attention imme- 
diately. The poor condition of much of the area is a result of past 
abuse by livestock and wildlife, compounded by steep topography, erosive 
soils, and occasional severe storms. At present, proper grazing manage- 
ment is probably the key to range improvement in the area. However, 
recreational and other uses may cause different abuses and problems which 
must be controlled. 


The major objectives of the research on the Reynolds Creek Watershed 

are to provide meaningful data for guidance in planning and implementing 
a range management and improvement program and to formulate procedures 
and models for predicting the effects of management and land-use changes 
on forage production, water yield, soil erosion, and related resources. 


During the past year, rest-rotation grazing management plans were drafted 
and discussed by BLM personnel, ranchers, and the Watershed Staff; however, 
the ARS research plans on the Reynolds Creek Watershed were not changed. 
Logging on private land was closely observed for evidence of soil erosion, 
but erosive storm events did not occur. 


Data collection for the 5-year study on the Rabbit Creek Watersheds was 
completed September 30, 1973, and instruments were removed soon thereafter. 
The final report on this study will be submitted about June 30, 1974, 

after analyses and data interpretations are complete. 


Data collection, processing, and analysis continued according to the 

F.Y. 1974 work plan and details of progress and accomplishments are G 
described in each section of the report. Further information is con- 

tained in Northwest Watershed Research Center Annual Reports for 1972 

and prior years and in Interim Report Nos. 1, 2, and 3 of ARS-BLM 

studies in the Reynolds Creek Watershed under Cooperative Agreement 

No. 14-11-0001-4162(N). 
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PRECIPITATION 


Title: Precipitation characteristics of a northern, mountainous 
semiarid watershed. 


Personnel Involved: 


Joes LUZ@l. Hyco les Lect. Supervise data collections and 
conduct special analyses requir- 
ing the use of a computer. 


W. J. Rawls, Hydrologist Supervise data reduction and 
: conduct special analyses requir- 
ing the use of a computer. 


Date of Initiation: June 1961 


Expected Termination Date: Continuing 


INTRODUCTION 


No dense, recording rain gage network existed in the Northwest prior to 
establishment of the Northwest Watershed Research Center. Such a net- 
work is necessary to delineate thunderstorms and storm variability. 
National Weather Service data collection stations are generally located 
in or near the main cities. Since these are generally along the main 
stems of major streams or in valleys, a sampling of precipitation on 
the range watershed areas is not available from their records. Also, 
there are too few rain gages capable of recording intensities or even 
individual storm data. 


Objectives: 


1. To develop methods for evaluating precipitation rates and amounts 
for watersheds of different sizes. 


2. To determine seasonal distribution of precipitation with respect 
to amounts, character, and areal extent. 


3. To develop depth-duration-frequency and depth-area-duration relation- 
ships for the Reynolds Creek Experimental Watershed. 


4. To establish general precipitation-elevation-aspect-slope- relation- 
ships from precipitation data obtained in the Reynolds Creek 
Experimental Watershed. 
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PROGRESS 


Precipitation data collection was continued from the gage network in 

the Reynolds Creek Watershed, comprised of 49 dual-gage sites (98 gages), 
and from the Rabbit Creek network, comprised of 5 dual-gage sites (10 
gages). The Rabbit Creek network was dismantled in October 1973. Fre- 
quency analysis and stochastic models were investigated as tools for 
describing temporal and spatial variations of precipitation. Using the 
dual-gage model, computer programs were developed to convert shielded 
and unshielded precipitation data to actual precipitation. In addition, 
a precipitation intensity program has been developed, and all precipi- 
tation data for 1968-1973 were converted to actual precipitation. The 
sensitivity of the dual-gage model coefficient was tested over its 
expected range and found to be relatively insensitive, producing error 
ranges of 10 percent or less. A comparison of rain gage data from gages 
equipped with free-swinging and constrained windshields indicates that 
gage catches are essentially the same for both types. 


The results of a comparison of annual unshielded gage catch and actual 
precipitation computed from the dual-gage model and regression equations _ 
are shown in Figure 1. The need for correcting unshielded gage catches 
is immediately apparent. For example, the unshielded gage catches of 
water years 1972 and 1973 are almost identical, yet a difference of 33 
(13 inches) cm of water exists in the computed precipitation values. 

This difference is attributable to less severe winds in 1973 to which 

the dual gage model is sensitive. 


A comparison of annual precipitation, computed from the dual-gage model 
and from a regression model, is shown in Figure 2. Regression model 
input consists of shielded gage values only, while dual-gage model input 
requires both shielded and unshielded gage values. The magnitude of the 
differences in precipitation suggests that a regression equation applied 
to either a free-swinging or constrained shield rain gage can produce 
precipitation values as reliable as the dual-gage model on an annual 

and seasonal basis. However, development of the simplified procedure 
was dependent on data from the dual-gage network. 


SIGNIFICANT FINDINGS 


The dual-gage precipitation model coefficient, when varied over its 
expected range, has proven to be relatively insensitive and produces 
error ranges of 10 percent or less. 


No difference in rain gage catch was found between rain gages equipped 
with constrained windshields and rain gages equipped with free-swinging 
windshields. 


A correction factor, derived from regression analysis'and applied to 
either a free-swinging or constrained shield rain gage, can produce 
precipitation values as reliable as the dual-gage precipitation model 
on an annual and seasonal basis. Since water balance and hydrologic 
model studies require an accurate precipitation input, computed actual 
precipitation can now more confidently be used as input to these models. 


WORK PLAN FOR FY 1975 


Frequency analysis and stochastic models will be used to determine and 
simulate temporal variations of precipitation. Studies of precipitation 
intensities will be initiated. 


REPORTS AND PUBLICATIONS 
PAUZeLs eee ee 
ABSTRACT - Estimation of rain gage catch deficiency. Pub. Trans. 
NW Reg. Meeting, Amer. Geophys. Union, Missoula, MT055.(2) 233s 
Pe Spero res 1973 3 


Precipitation studies at Reynolds Creek. Pres. ARS-SCS Western 
Hydrol. Workshop. Boise, ID. 
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Unshielded gage precipitation catch and computed actual precipitation, 
Station 176X07, Elevation 6600 ft. 
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SNOW 
Title: Snow accumulation, snow redistribution, and snowmelt. 
Personnel Involved: 

L. M. Cox, Hydrologist Supervise the planning, designing, 
execution, analyzing, and report- 
ing of proposed experiments. 

J. F. Zuzel, Hydrol. Tech. Assist in the planning, designing, 
execution, analyzing, and report- 
ing of proposed experiments. 

W. J. Rawls, Hydrologist Assist in the planning, designing, 


execution, analyzing, and report- 
ing of proposed experiments. 


Daterottinit tations 72.964 


Expected Termination Date: Continuing 


INTRODUCTION 


A substantial proportion of the runoff from the western rangelands has 
its origin in rapid melting snow. To improve the quantity or timing of 
flow from snow-fed streams by manipulation of vegetation or by other 
practices requires that the behavior of snow be well understood. There 
has been littlé research on the behavior of snow in shrub areas any- 
where-——and almost none in the sagebrush areas of the Northwest. 


Destructive late winter and spring floods in the Northwest frequently 
originate from rapid melting snow at low elevations, characteristic of 
the sagebrush zone. Although there is little likelihood of modifying 


‘snowmelt rates enough to alleviate this threat, knowledge about the 


behavior of snow in the sagebrush zone will be helpful in devising 
better warning and forecasting techniques that may reduce the danger 
to life and property from snowmelt floods. 


Objectives: 


1. To determine the physical and meteorological factors contributing 
to nonuniformity of snow accumulation in shrub-covered study 
basins on mountainous terrain. 


2. To determine the influence of the controlling physical and meteor- 
ological factors on snowmelt from the above areas. 


3. To improve snowmelt prediction techniques by evaluating the energy 
exchange process of the snow surface under different snow cover 
conditions. 
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to minimize evaposublimation losses and to maximize and prolong 


4. To study the oasis effect of isolated, late-lying snowdrifts € 
water yield from snowmelt. 


PROGRESS 
Snowmelt evaluation studies showed that vapor pressure, net radiation, 
and wind explained 78 percent of the variance associated with the 
snowmelt process. Melt predictions, based on vapor pressure, net 
radiation, and wind, were improved by 14 percent as compared to using 
temperature measurements alone. 


A long-term water supply forecast model, using parameter optimization 
techniques on snow course data, was developed and successfully tested 
on three basins in southwest Idaho. 


Preliminary analysis of photogrammetric snow measurements indicates 
that snow depths in the study area can be determined with a standard 
error of 0.6 ft. Snow density zoning methods were not critical in 
estimating the total quantity of water stored in a small watershed 
snowpack. 


SIGNIFICANT FINDINGS 





Snowmelt predictions, based on vapor pressure, net radiation, and 
wind, improved melt predictions by 14 percent, compared with predictions 
using temperature measurements alone. 


A long-term water supply forecast model has improved forecast accuracy 
by 5 percent when compared to the models now used by the Snow Survey 
Section of SCS. 


WORK PLAN FOR FY 1975 


The energy exchange process will be evaluated by comparing data from 
meteorological parameters, independent measurements of energy balance 
components, and melt collectors. The data were obtained during FY 1974 
on areas of continuous snow cover and isolated late-season snowdrifts. 


Short-term streamflow forecasting techniques will be further refined, 
using melt data from several types of snowmelt collectors, snow courses, 
and energy balance measurements. 


The accuracy of predictions of the areal distribution of snow and snow- 
water content on a watershed will be tested, and 3 years of photogrammetric 
and 10 years of snow course data will be statistically evaluated. 


bef 


REPORTS AND PUBLICATIONS 


Costa. Me i975 
ABSTRACT - A new device for evaluating the energy exchange 
processes that contribute to the ablation of a snowpack. 
20th Ann. Meet. Pacific NW Reg. Sec., Amer. Geophys. Union. 55(2):77. 


COX nce eee se uzel L975 
Forecasting runoff from universal surface gage snowmelt measure- 
ments. Jour. of Soil and Water Cons. 28(3): 131-134 


Norris, vo. eo; .eandyLe M. “Cox 1974 
Research plus action. Soil Cons. 39(6). 


COx, ise, eis 
Snow hydrology studies at the Northwest Watershed Research 
Center. Pres. ARS-SCS West. Hydrol. Workshop, Aug. 20-24, 
Boise, ID. 


Rawls, We J.,,andd. R. Zuzel 1973 
Computer modeling studies. Pres. ARS-SCS West. Hydrol. Workshop 
Aug. 20-24, Boise, ID. 


Cox, oer. 
A device for evaluating the net water vapor exchange between 
snow and air. (Being reviewed for pub. in Water Resources Res.) 


ZU2Z@) 0 J3.0F., anid Ls. My, Cox 
A statistical analysis of the relative importance of meteoro- 
logical variables in the snowmelt process. (Being reviewed for 
pub. in Water Resources Res.) 


Zuzel, J. F., D. C. Robertson, and W. J. Rawls 
Optimization of long-term streamflow forecast model coefficients. 
(Being reviewed for pub. in Jour. of Soil and Water Cons.) 
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VEGETATION AND SOIL MOISTURE 


Title: Evaluation of cover production, herbage yield, and soil 
conditions for different levels of vegetation management. 


Personnel Involved: 


G. A. Schumaker, Soil Scientist Plan, design, and coordinate 
research activities and prepare 
reports. 

W. J. Rawls, Hydrologist Prepare computer programs for 


analysis of vegetation data. 


C. L. Hanson, Agr. Engin. Perform computer analysis relative 
to soil moisture data and assist 
in summarizing field data. 


Digi GOOD, v¥orol. Lack. Responsible for various aspects of 
data collection and field obser- 
vations including soil moisture 
Measurement and calibration; 
compile research outline data. 


Date of Initiation: May 1971 


Expected Termination Date: Continuing 


INTRODUCTION 


Quantitative data on herbage yield from rangelands under different levels 
of management are needed to guide land managers in coordinated multiple 
use of the range. These needs require more discerning information on how 
vegetation and soils respond to imposed treatments, including controlled 
grazing. Information is also needed with regard to methods of increasing 
cover and to the rate of recovery of native range following intensive 
grazing practices. 


Objectives: 


1. To determine the effects of grazing management and treatments on 
yield of herbage, cover production, soil moisture regimes, and 
soil surface conditions at selected sites. 


2. To study changes in plant density and plant composition as a result 
of grazing management and treatments. 


é 
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PROGRESS 


Vegetation 


Vegetation studies for 1973 were concentrated on six of the nine study 
sites at the Reynolds Creek Experimental Watershed. Cover changes occur 
slowly on the other three sites that have sparse vegetation; therefore, 
clippings and basal area data were not obtained. Descriptive infor- 
mation for all study sites is given in Table 1. Locations of the study 
sites are shown on Figure 1 of the Introduction section of this report. 


Herbage yields: The double sampling technique was used once again this 
year in determining plot yields. This method involved two steps: 

(1) three sample areas, which previously had been estimated, were clipped 
and weighed in the field by two technicians; (2) herbage weights were 

then estimated on ten other random sample areas. This procedure increased 
the sample size with only a minimum of extra time. 


Total herbage and non-sage yields from the brush treatments at the Nancy 
site, where the annual precipitation is about 11 inches, are shown on 


Figure 1. Some small sagebrush plants were evident on the brush-removed 


treatment, which accounts for the 200 lb/ac difference between the total 
herbage and non-sage yields. These plants were in the seedling stage at 
the time of clearing in 1971. The herbaceous sagebrush on the sprayed 
treatment contributed about 150 lb/ac, where about 80 percent kill was 
accomplished by spraying in 1971. A statistical analysis showed that 
the yield of grasses and forbs on the brush-controlled treatments was 
significantly greater than on the other treatments at the 5 percent level. 
Competition for moisture and shading by the brush limited the yields of 
grass and forbs where the brush was controlled. It is speculated that 
a new stand of grass was slow to establish at this site following spray- 
ing in 1971. 


Yields in 1973 from the Upper Sheep Creek study are shown on Figure 2. 
Total herbage yields on all three treatments within the exclosure were 


. near 2,000 lb/ac. There was a very small amount of Artemesia tridentata 


on the sprayed and brush-removed areas after five growing seasons 
following treatment. About 1,000 1lb/ac more grasses and forbs were 
produced on the sprayed and the brush-removed treatments than were 


produced without brush treatment. This difference was statistically 


significant at the 5 percent level. 


Total yields and non-sage yields at the Reynolds Mountain site are 
shown on Figure 3. Grass and forb yields on the brush removed area 
were significantly greater at the 5 percent level than on the other 
treatments. It is apparent that more than 2 years response time is 
needed for the understory grasses and forbs to increase following spray 
treatment. Sagebrush on this site was sprayed in 1971. 
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The Whiskey Hill site includes a spray treatment but not a brush-removed 
treatment. Total yields in 1973 were 503, 836, and 800 lb/ac on the 
sprayed, untreated, and grazed treatments, respectively. Only one 
growing season has passed since the sagebrush was sprayed, and total 
yields on the sprayed treatment have not yet equaled or surpassed the 
yields on either the grazed or untreated sites. Non-brush yields were 
493, 449, and 302 1b/ac on the sprayed, untreated, and grazed treatments, 
respectively. More grasses and non-brush forbs were produced on the 
sprayed treatment than on either of the other two treatments. The in- 
crease was statistically significant at the 10 percent level. 


The Nettleton site has an excellent cover of grasses and forbs with only 
a small percentage of brush species present. Before 1971 the site was 
under private management and was grazed late in the season. Since the 
site was developed for experimental purposes, the grazed portion of the 
study has been subjected to intensive grazing. Total yields in 1973 
were 1,133 lb/ac from the exclosure and 1,208 lb/ac from the grazed 
portions of the study. Although vegetation on the Nettleton study site 
appears uniform, yields from individual sample points ranged from 500 to 
1,500 1lb/ac in the exclosure. 


The Nettleton site was subjected to intensive grazing by 12 cattle during 
the 21-day period from June 11 to 25, 1973. The 6.33 acres provided .53 
acres per animal for grazing. From the clippings, the forage production 
was estimated to be 1200 lb/ac and should have been nearly sufficient 

for animal requirements during the grazing period. 





Species composition: Vegetation data from the 700 line transect points, 
taken on each of the brush treatments at the four brush treatment sites, 
were divided into four groups: grasses, lupine, other forbs, and shrubs. 
The study at Upper Sheep Creek was started in 1969. The 1973 information 
for this site is of interest for comparison of different treatments. The 
percent of various vegetation species groups is shown on Figure 4; hori- 
zontal lines provide a graphic comparison of response to different treat- 
ments. The small percentages of brush species on either the brush-removed 
or sprayed treatments show the effectiveness of brush control. The percent 
of lupine was not reduced by spray application, but other forbs were 
controlled by spraying. ; 


Comparisons of similar groups of vegetation at the Reynolds Mountain 
site are shown on Figure 5. Lupine made up a large portion of the 
vegetation on the brush-removed treatment. This was probably due to 
the mechanical nature of the grubbing operation, which resulted in 
better seed placement. At present, groups of vegetation are not 
different on either the grazed area or the untreated portion of the 
exc losure. 
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Similar data are given in Table 2 for the Nancy and Whiskey Hill sites. 
No distinct trends are apparent at this time other than the effect of 
the treatments on brush species. 


Vegetative cover: Line transects for depicting cover were run in 1973 
on all but three of the study sites. The percentage for each cover 
class at each of the sites is given in Table 3. In general, the grazed 
treatments have less vegetative cover and more bare ground than the 
untreated portion of the exclosure. 


The sprayed treatment at Reynolds Mountain has a comparatively low 
amount of vegetation and more bare ground. Continued measurements will 
determine changes in cover makeup due to treatment. 


Trend plots: A series of photographs were obtained on the trend plots 

at all study sites to depict changes in vegetation under the different 
treatments. These photographs have been taken annually from the beginning 
of the study. 


Soils 


Erosion condition: Observations for soil Surface Factor determinations 
were taken on the different treatments at each of the study sites in 1973. 
Very little change was noted from 1972. All sites were given the erosion 
condition class of "stable" except a designation of “slight"°for the 
sparse vegetation sites at Upper Sheep Creek and Lower Sheep Creek. 


Comparison of bulk densities: Bulk density measurements were taken at 
Whiskey Hill following grazing to determine differences between the 
exclosure where cattle did not have access and the adjacent area open 

to grazing. Average bulk densities from ten sampling points were becv 
g/cc on the sprayed area of the exclosure and 1.30 g/cc on the grazed 
area. Statistical analysis showed that these means were significantly 
different at the 5 percent level. The difference is attributed mainly to 
the effect of cattle trampling. There appears to be a sufficient increase 
in grass to cause increased root permeation in the soil at the time of° 
measurement, which partly accounts for the decreased bulk density. 


Soil water: Soil water measurements were taken regularly during the 
growing season, beginning when the soil water was near its maximum and 
ending when it was depleted. In comparison with the other treatments, 
brush removal showed some effect on water loss at the three sites: Nancy, 
Upper Sheep Creek, and Reynolds Mountain (Figures 6.av,nand,6). There 
was very little difference in water loss between the sprayed and untreated 
plots. The cover offered by the dead sagebrush on the sprayed treatment 
probably acted much the same as the live sagebrush in providing shade to 
the soil surface and reducing air movement, thus keeping surface soil 
moisture loss to a minimum. 


At Reynolds Mountain a comparison was made of water removed on the 
three different brush treatments from the 4-foot zone to determine ¢ 
water use by sagebrush late in the season. There were no apparent 
differences in water use between treatments. August and September 

rain at the Reynolds site may have provided necessary moisture for 

growth rather than the sagebrush drawing moisture from the 4-foot 

zone. 


Daily evapotranspiration rates were determined for the period of June 

5 through July 31 at Upper Sheep Creek, including precipitation during 
this period. Similar data was summarized for an 84-day period 
beginning May 9 at Reynolds Mountain. At both locations the daily 

loss due to evapotranspiration was .09 inches of water per day. This 
value agrees very favorably with values reported by Rawls (see publi- 
cation list at end of this section) where estimates from data collected 
over a period of years showed .10 in/day as evapotranspiration at 
Reynolds Mountain and .07 in/day at Lower Sheep Creek. 


SIGNIFICANT FINDINGS 


At Upper Sheep Creek and Reynolds Mountain, the amount of non-sage 
herbage was significantly greater on the sprayed and brush-removed 
treatments than where the sagebrush had been treated. 


The percent of lupine was much greater where sagebrush was removed by f 
grubbing in 1971 at Reynolds Mountain. 


Significant differences in soil bulk densities, measured on grazed and 
ungrazed sites at Whiskey Hill, were attributed mainly to the effect of 
cattle traffic during grazing. 


Moisture loss from the soil was different where Artemesia tridentata 
remained than where the brush was completely removed, leaving only 
understory grasses and forbs. 


WORK PLAN FOR FY 1975 


Measurement of changes in species composition, cover, herbage yield, 
and soil conditions for different vegetative treatments will be 
continued. Emphasis will be given to measurement of herbage yield, 
notation of cover change, and measurement of changes in plant vigor 
at the different sites. 


Yield and vegetative cover data collected in the field will be related 
to treatments on sites representing different ecological zones. 
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« Water use under different types of brush control will be investigated 
») by measurement of neutron soil moisture through the growing season at 
five study sites. 


Methods will be developed for sampling the surface inches of soil on 
sites where gravelly soils do not lend themselves to ordinary methods 
of bulk density measurement. 


The measurement of change in crown height of sod due to frost heaving 
will be continued at two sites. 


The percent of bare ground at the end of the grazing season will be 
measured to determine forage utilization and erosion potential. 


The effect of cattle traffic on soil compaction and the effect of 
increased grass cover on bulk density will be studied. 

REPORTS AND PUBLICATIONS 
Rawls, W. J., et al. 1973 


Soil moisture trends on sagebrush rangelands. Jour. of Soil and 
Water Conserv. 28(6): 270-272. 
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TABLE 1,--Tabulation of site information 











2/ ad EGETATIVE2/ 
ELEVATION| SLOPE ASPECT PRECIPITATION 2/ COVER | SCS HYDROLOGIC 
SITE (feet) (percent) OF SLOPE (inches) SOIL SERIES ~ | (percent)] CLASSIFICATION 


SPARSE VEGETALIN SITES 


Flats 4000 Be North Nannyton loam 
Upper Sheep Creek 6100 Southwest pravelly loam 
Reynolds Mountain ater Bullrey 

West 6850 


Southwest jae pravelly loam <2 - &B 
Nettleton's _ 5000 
te Bullrey 
32% pravelly loam B 
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1/ Good grass cover site on soil developed from basalt. 
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2/ Based on 1961 survey. 
3/ Snow removed by wind. 


4/ Snow deposition zone. 
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TABLE 2.--Vegetative composition at Whiskey Hi112/ and Nancy study } 
sites expressed as percent of total vegetation. ¢ 


Whiskey Hill 


Vegetation 


Grasses 
Lupine 
Other Forbs 
Shrubs= 





Nancy 


Vegetation 


Grasses 
Lupine 


Other Forbs 
Shrubs2/ 





1/ There was no brush removal treatment at the Whiskey Hill site. 


2/ Including Artemesia tridentata 
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- TABLE 3.--1973 cover measurements taken at study sites. 


Brush Removed 
Sprayed 
Untreated 
Grazed 
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Brush Removed 
Sprayed 
Untreated 
Grazed 


Upper Sheep Creek 





Brush Removed 
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Figure 7. Soil water content for the O- to 3-ft. depth during the growing season and 
daily precipitation at the Upper Sheep Creek site for three different brush 


treatments. 
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INFILTRATION 


Title: Developing, testing, and evaluating an analytical infiltration 
model 


Personnel Involved: 


W. J. Rawls ~ Supervise field data collection; 

(ARS-Boise) coordinate the research with the 
cooperators; perform any other 
analysis needed. 


G. A. Schumaker Collect necessary soils and vege- 
(ARS—Boise) tation data 

W. R. Hamon Coordinate the 5-year technical 
(ARS-Coshocton, Ohio) bulletin on infiltration work per- 


formed at Boise 


R. W. Jeppson Develop and field test analytical 
(Utah State University) infiltration models 

R. H. Brooks Assess soil properties in laboratory 
(Oregon State University) and in the field; help establish 


field procedures; help develop and 
field test analytical infiltration 
models 


Date of Initiation: February 1971 


Expected Termination Date: Continuing 


INTRODUCTION 


Infiltration of water into soil profiles is probably the most critical 
hydrologic component in watershed management, overland flow prediction, 
sediment generation, natural and artificial ground-water recharge, and 
irrigation. More accurate information concerning the flow system result- 
ing from infiltration is essential if agricultural lands are to be managed 
for optimum multiple use. The Northwest Watershed Research Center has 
been conducting an infiltration research program aimed at developing, 
testing, and evaluating analytical infiltration models. 


Objectives: 


1. To develop or adapt mathematical models in the form of partial 
differential equations to describe steady-state and transient 
one-dimensional and three-dimensional axisymmetric flow through 
partially saturated soils. 


4 =e. gates ’ 
et SSS a = weer em ee ee ae Dees rele a One es, BAe g Bt eee see Me Seem nn wt ee a en Oe 


Pap To use mathematical models to adjust for the lateral "spreading € } 
effect" of moisture movement from a circular rainfall simulator 
under various soil types and conditions. 


a. To test and refine the mathematical models by comparing results 
with laboratory and field determinations of infiltration, and to 
establish the relative influence of the several interacting 
physical processes on infiltration. - 


4. To determine the saturation-pressure relationship (moisture 
characteristic) and saturated conductivity parameters by the use 
of parameter optimization in conjunction with the mathematical 
models and infiltration data. 


ae To determine the potential quantity of water retainable by various 
soil-vegetation complexes independent of infiltration and different 
initial soil moisture levels. 


Mathematical models of steady and transient one-dimensional and three- 
dimensional axisymmetric flow through partially saturated soils have 
been formulated by R. W. Jeppson of Utah State to satisfy objectives 
No. Ll and 2. The collection of the field data to satisfy objectives 
No. 3, 4, and 5 was performed in the summer of 1972. The field pro- 
gram was discussed in the Interim Report No. 3. 


PROGRESS 


The numerical infiltration model testing with the field data was 

expanded beyond the models developed by Jeppson at Utah State University, 
to include those proposed by Rogowski (ARS-Pennsylvania State University), 
Brooks (Oregon State University), and Smith (ARS-Tucson, Arizona). The 
results from this expanded testing program have not been completed; 
however, preliminary results indicate that all the models yield satis- 
factory results. Testing of a number of the well known empirical 
infiltration models has also begun. After the testing of all the models 
has been completed, they will be evaluated according to their useful- 
ness in evaluating the effects of land management changes. 


Laboratory capillary pressure-saturation analysis of soil cores taken 
during the field program were completed. Comparison of insitu capillary 
pressure-saturation data and laboratory capillary pressure-saturation 
data indicated a good correspondence at the surface or near the surface 
(0 to 2 inches); however, at lower depths the correspondence was very 
poor. This poor correspondence was due to sampling procedures and 
sampling conditions. 
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The numerical models indicate that the surface layer has the most 
predominant effect upon the hydraulic behavior of the soil. Thus, 
hydraulic classification of the surface soils may be all that is 
necessary to describe the infiltration characteristics of the soil. 


Attempts at fitting a number of empirical equations describing the 
capillary pressure-saturation relationships under drainage and 
imbibition conditions indicate that the imbibition and drainage 
functions are different. Most of the present equations were developed 
for the drainage conditions. Our research indicated a more suitable 
function for imbibition was: 


Sa = -85b 
(Po/P,* ay 
where, 
S. = effective saturation 


P = capillary pressure 


bubbling pressure 


rd 
o 
i] 


X = pore size distribution 


b constant 


The unknowns in this function are b, PS and A and attempts are being 
made relating these to soil properties. 


A procedure for determining the hydraulic properties of natural soils 
insitu was developed and can be summarized as follows: 


1. Soil density and capillary pressure data were obtained simultan- 
eously at designated depths with an infiltrometer, gamma probe, 
and tensiometer equipment through the time that water was 
applied at a constant rate less than the ultimate infiltration 
capacity of the soil. 


2. Bulk densities of soil samples taken from the site were deter- 
mined to convert the soil density measurements to saturation 
values.. 


3. Numerical solutions simulating the field tests and based on the 
field saturation-capillary pressure data from steps 1 and 2 
were obtained to define the unsaturated hydraulic properties of 
the soil. 
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4. The optimal dimensionless application rate was determined by ¢ 
matching the numerical solutions to the saturation-time curve 
of the field data. The saturated hydraulic conductivity was 
then computed. 


oi Upon achieving satisfactory agreement in Step No. 4, the 
mathematical model could be used in sensitivity analyses to_ 
study such effects as different application rates, initial 
conditions, etc., on moisture spreading, rate of movement 
and location of wetting front, and other factors associated 
with the infiltration process. 


Four chapters of the technical report on the infiltration work at Boise 
were completed. Work is continuing on the last two chapters. 


SIGNIFICANT FINDINGS 


A procedure for determining the hydraulic properties of soils insitu 
was developed. It was found that present equations representing the 
capillary pressure-saturation relationship during drainage are not 
applicable for imbibition. Also, a comparison of insitu capillary 
pressure-saturation data:and laboratory capillary pressure-saturation 
data indicated a good correspondence at the surface or just below the 
surface; however, this was not true at lower depths. 


WORK PLAN FOR F.Y. 1975 | 


The relative influence of the physical processes on infiltration will 
be studied using optimization, and various theoretical and empirical 
infiltration models will be field verified. A technical bulletin 
summarizing all the infiltration work at Boise will be completed. 


REPORTS AND PUBLICATIONS 


Brooks, Royal H., Paul J. Leclercq, Richard R. Tebbs, and Walter J. 
Rawls 1974 
Axisymmetric infiltration. Water Resources Res. Inst. Report 22, 
Oregon State Univ., Corvallis, OR. 61 pp. 


Rawls, Wi Jeo 1973 
Infiltration studies pres. at ARS-SCS Western Hydrol. Workshop, 
August 20-24, Boise, ID. 


Will, Gregory R. 1973 
Rainfall simulator. Masters Thesis, Oregon State Univ., 70 pp. 
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EVAPOTRANSPIRATION 


Title: Natural evaporation from sagebrush rangelands, alfalfa, and 
stock ponds in a semiarid environment. 


Personnel Involved: 


C. L. Hanson, Agr. Engin. Plan programs and procedures; 
; design and construct facilities 
for evaporation studies. Perform 
analyses and summarize results. — 


G. H. Belt, Assoc. Prof. Analyze and report on experimental 
Watershed Mangt. data. 
University of Idaho 


Date of Initiation: November 1968 


Expected Termination Date: Continuing 


INTRODUCTION 


A complete understanding of the evaporative process is essential when 
developing predictive relationships for the evapotranspiration com- 
ponent in definable soil vegetation complexes under a particular level 
of management. The Northwest Watershed Research Center, in cooperation 
with the Department of Forestry, University of Idaho, has been conduct- 
ing evapotranspiration studies designed for measuring and predicting 
evapotranspiration under sparse vegetative cover and unsaturated 
surfaces. 


Objectives: 


1. To determine the evaporative loss of water from sagebrush range- 
lands, irrigated alfalfa, and stock ponds, while observing 
pertinent meteorological parameters and the soil moisture status. 


2. To develop, for predictive purposes, relationships for associ- 
ating the evaporative loss with meteorological parameters, type 
and degree of surface cover, soil moisture, and potential 
evaporative demand. 
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PROGRESS 


Profile micrometeorological data obtained over irrigated alfalfa and 
dryland sites with sparse covers of sagebrush and greasewood were 
utilized to determine the Bowen ratio and the surface and aerodynamic 
resistances. The surface resistance was found to be proportional to 
the Bowen ratio. The proportionality coefficient, a bulk surface 
resistance parameter, was dependent upon the availability of water for 
the evaporation process. This resistance representation of the Bowen 
ratio was combined with the aerodynamic and energy budget equations 

to obtain an energy balance-resistance model; an iterative procedure 
was used to compute evapotranspiration and near-surface temperature. 
The data required consists of air temperature and humidity, windspeed, 
surface roughness, and an estimate of the bulk surface resistance 
parameter. This parameter was incorporated into a combination equation, 
and computed values of evapotranspiration were compared with those 
obtained by the energy balance-resistance model. Development of this 
model is continuing at the University of Idaho. 


A somewhat different approach from the previous work is now being used 
to describe watershed evapotranspiration. The first reason for this 
change in direction is because of the change in personnel; and the 
second being an attempt to develop an evapotranspiration model that 
can be directly included in a watershed hydrologic model. 


A review of the literature indicates that a model of the following 
form can be developed for the Reynolds Creek Watershed. 


The model is expressed as: 
ET = ET, KK sw + EV 
in which 


ET 


total evapotranspiration (inches/day) 


eo 
tr 
il 


2 potential evapotranspiration (inches/day) 


A 
i] 


plant coefficient 


coefficient of limiting soil-water 


A 
) 


EV = evaporation from the soil and plants (inches/day) 


Four lysimeters, two at the Lower Sheep Creek and two at the Reynolds 
Mountain experimental sites, were calibrated and the necessary instru- 
ments obtained so that these lysimeters will be in operation during 
the 1974 growing season. The information from these lysimeters will 
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be used to obtain the constants for the above evapotranspiration model. ¢ 
Data from the neutron soil-water measuring program will also be used 

to improve the model and to determine if the constants in the evapo- 
transpiration model will have to be different for various areas on the 
watershed. 


| SIGNIFICANT FINDINGS 
An energy balance-resistance model for computing evapotranspiration 
was developed by using profile micrometeorological data obtained over 
irrigated alfalfa and sagebrush rangeland. Data required for the model 
consists of air temperature and humidity, windspeed, surface roughness, 
and an estimate of the effective resistance of the surface. 


WORK PLAN FOR FY 1975 


The goal is to develop a basic evapotranspiration model for rangeland 
watersheds, based on the equation 


et Pek hr Bs 
Pp p sw 

using lysimeter data and supporting information from the neutron soil- 
water measuring program. This model has been successfully used to 
describe evapotranspiration from rangeland in South Dakota. It will é 
most likely be necessary to obtain additional information from the , 
lysimeters so that the constants in the model can be tested. The 
model will then be used in a watershed hydrologic model. 


A manuscript describing the bulk energy balance-resistance evapo- 
transpiration model will be completed. 


REPORTS AND PUBLICATIONS 


Hamon, W. R., and G. H. Belt. ABSTRACT. Energy balance-resistance 
model for computing evapotranspiration. Trans., American Geophys. 
Union, 54(4) H59: 270, 1973. 


Hanson, C. L. Model for predicting evapotranspiration from native = 
rangelands in the northern Great Plains. Ph.D Dissertation, Utah 
State Univ., Logan, UT, 116 pp, 1973. 
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WATER QUALITY 


‘Title: Water quality characteristics of the hydrologic flow regime 
of the Reynolds Creek Experimental Watershed. 


Personnel Involved: 


G. R. Stephenson, Geologist Responsible for coordinating 
activities with cooperators. 
Design collection network and 
responsible for project completion. 


J. F. Zuzel, Hydrol. Tech. Responsible for statistical analysis 
of data and shares the responsibility 
for aquatic sampling. 


C. M. Rountree, Biol. Tech. Responsible for collection of water 
samples, aquatic samples, and labora- 
tory analyses. 


‘Date of Initiation: October 1972 


Expected Termination Date: Continuing 


INTRODUCTION 


In recent years, because of the increased concern for the quality of our 
environment, many agricultural practices have come under close scrutiny 
as potential sources of air and water pollution. 


Several agricultural practices are known to have contributed to pollution 
of surface and ground-water resources. Feedlot operations under certain 
conditions do contribute to increasing nitrate content of adjacent surface 
and ground water. Heavy grazing on open range usually compacts the soil, 
reduces infiltration, and seriously reduces the vegetative cover, increas- 
ing runoff. Increased turbidity and contribution of sediment and nutrients 
to the streams are often the result. Irrigation usually contributes 
heavily to total salt and nitrate-nitrogen content downstream. However, 
not all agricultural operations are detrimental. 


Information is needed on the water quality characteristics of rangeland 
watersheds under natural conditions and various land management practices. 
The Reynolds Creek Experimental Watershed offers an excellent opportunity 
to study water quality characteristics related to several of the above- 
mentioned problems. No commercial fertilizers or pesticides have been 
used on the watershed. Herbicides were used infrequently for sagebrush 
control, but not since 1965. 


; - ‘ ; 5 si 
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. With the present distribution of hydrologic networks throughout the ¢ 


watershed, sampling of both surface and subsurface flow for water 
quality analyses can easily be accomplished. The water quality con- 
gtituents can be related to the hydrology of the system, particularly 
the properties of the water, the distribution, and the circulation. 


The BLM has expressed need for more information on water quality changes 
influenced by various land management practices. As more rangeland is 
being used for recreation, this information becomes more important. 


Objectives: 


To determine water quality characteristics of the hydrologic flow 
regime of the Reynolds Creek Experimental Watershed as related to 


1. Concentrations of cattle on local areas of rangeland and 
quasi-feedlot conditions, 


2. Irrigation return flow, and 


3. Natural soil and geologic conditions. 


PROGRESS 


During the past year, seven samplings for water quality have been made & 
at nine sites on Reynolds Creek. Approximately 25 chemical analyses 

were completed for each sample, providing baseline data for determining 

water quality characteristics. The samples were also analyzed for total 

and fecal coliform counts. 


Four samplings of aquatic invertebrates were made during the year at 
five sites on Reynolds Creek, ranging in elevations from 3600 to 5400 
feet. Three samples were taken at each of the five sites using a 
Surber square-foot sampler. Samples were preserved in alcohol, sorted 
to order, and counted. 


Because of the near drought conditions existing on Reynolds Creek 
throughout most of 1973, any conclusions reached from analysis of this 
data should be related only to unusually low flow conditions. The 
prospects for greater streamflow in 1975 are much improved. This 
should aid in observing quality variation under different climatic 
conditions. 


Chemical analysis: Because of the quantity of data collected during 
the water year at all nine sites on Reynolds Creek, only two selected 
sites will be used here to illustrate the significant water quality 
variations that occurred. The Lower Reynolds site characterizes waters 
receiving irrigation return flows, and the Tollgate site characterizes 
streamflow from open rangeland, upstream from any irrigation return flow. e 
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In comparing the dominant chemical constituents at the two sites 
(Figure 1), there is a marked increase in levels of concentration at 
the Lower Reynolds site during the irrigation season. Concentration 
levels at the Tollgate site remained relatively constant. The sodium 
adsorption ratio (SAR), an alkali hazard index computed using sodiun, 
calcium, and magnesium concentrations, shows that the water at Lower 
Reynolds is not dominated by sodium and has a low alkali hazard, but 
that there is an increase in salinity during the irrigation season. 
This same relationship is evident when analyzing the electrical con- 
ductivity (EC) measurements. Other chemical constituents shown in 
Figure 1, though present in much lower concentrations, show much the 
same increase during the irrigation season. 


Total phosphate (P), often used as an indicator of eutrophic conditions 
because of the solubility of the orthophosphates from applied fertilizers, 
does not follow the same variation trends in Figure 1, as do the other 
constituents. Phosphate ions characteristically are adsorbed to the 
suspended solids in streamflow and to channel sediments. This condition 
appears to be true in Reynolds Creek because the higher the percent of 
suspended solids the higher the concentration of total (P). Because 

no commercial fertilizers have been applied to range or croplands on 

the Reynolds Creek Watershed, the total phosphates are derived from 
natural sources and very little soluble orthophosphates are present. 


From the above analyses, based upon a relatively dry year, the increase 
in chemical concentration progressing downstream on Reynolds Creek is 
caused by irrigation and very little change is related to grazing 
conditions on open range. 


Coliform analysis: Coliform analyses of the water samples were not 
conclusive because counts vary considerably between sites and between 
sampling dates. However, it appears that average fecal coliform counts 
are higher at sites where there are concentrations of cattle (See Table 1). 
Further analyses of the coliform data show that an increase in both 

total and fecal counts occurs with a corresponding increase in flow. 

A need for an increase in the sampling network and sampling frequencies 

is evident. ; 


TABLE 1.—Average fecal coliform counts for the Reynolds Creek 
sampling sites. 


Reynolds Outlet 78 | *Nettleton 166 Democrat 30 
*Lower Reynolds 104 Tollgate 56 |*Reynolds Mtn. Weir 181 


Tyson 89 Dobson Creek 86 Reynolds Mtn. Spring 8 





* These sites are subject to concentration of cattle. 
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Aquatic invertebrates: All aquatic samples were dominated by immature ¢ 
insecta with a scattering of leaches and snails, especially at the lower 
elevations. Total populations ranged from a high of over 1300 individ- 

uals per square foot at the outlet weir (3600 ft. msl) in October to 

90 per square foot at the Tollgate Weir (4600 ft. msl) in July. Average 
sample composition was: Trichoptera 47%, Ephemeroptera 23%, Diptera 154, 
Plecoptera 7%, and Coleoptera 6%. Invertebrate forms other than the 

insecta accounted for 2% of the average sample composition. 


Although the data are quite limited, tentative inferences as to the 
spatial and temporal variations in total population can be made at this 
point. The results of a temporal and spatial analysis of the data, 

using elevation as the space defining variable and month of the year 

as the time defining variable, are given on Figure 2. Time and elevation 
trends are evident. Generally, total populations are highest at all 
stations in October and lowest in June. Since August totals are higher 
than in June, additional sampling should define the peak of this apparent 
cyclical change. Channel relationships being equal, the elevation to 
population relationship indicates that less dense populations will be 
found at increasingly higher elevations. Quantification of this relation- 
ship will also require additional sampling. An examination of Figure Ze 
suggests the possibility that the combination of a log transform of the 
space parameter and an autoregressive model, using the time parameter, 
could define and possibly predict the population response surface. 


Ground water investigations: Geochemical investigations of multi-aquifer ¢ 
systems in basalt and lake sediments have provided data for selecting 

several chemical parameters for use in identifying the waters of two 
different aquifer systems. The relationship of sodium adsorption ratio 

to electrical conductivity was found to be the most significant. 

Differences in this relationship were used to classify the waters accord- 

ing to origin. Basaltic waters were found suitable in quality for most 

crops grown in this area; whereas, water from lake sediment aquifers was 

less suitable. Considerable variation in water quality exists where there 

is hydraulic continuity between the two aquifers. 


SIGNIFICANT FINDINGS 


Analysis of 1973 surface water data show that: (1) increases in chemical 
concentration of surface water in Reynolds Creek, such as increased 
salinity, are related to irrigation return flow and not to grazing con- 
ditions on open range, and (2) increases in total and fecal coliform 
counts are related to concentrations of cattle, but are still much 
below Public Health Standards. 


Geochemical investigations of ground-water flow systems show that: 
(1) several chemical parameters may be used to identify waters of 
separate systems, and (2) water from the basaltic aquifer is much 
better in quality for irrigation than lake sediment water. € 


. 2 
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WORK PLAN FOR FY 1975 


Additional sampling sites will be added to the present sampling network 
in areas where cattle concentrate during the grazing season on open 
range. These sites will be used to obtain better coliform information. 
Additional samples will be collected during high flow conditions when 
they occur. 


The present sampling network will be monitored for both chemical and 
aquatic determinations. 


We will continue to analyze the data as soon as possible after laboratory 
analyses are completed. No change in method of analyses, as outlined in 
the 1974 annual report, is contemplated. 

Water quality data from ground water and surface water systems will be 
used in a comprehensive water quality model. 


REPORTS AND PUBLICATIONS 


None. 
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Figure l. 
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Major chemical constituents (mgl) and channel flow at Tollgate 
and Lower Reynolds sites. 
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RUNOFF AND SEDIMENT 


Title: Sediment yield from rangeland watersheds. 


Personnel Involved: 





C. W. Johnson, Supvy. Hydraul. Plan programs and procedures; 

Engin. design and construct facilities 
for sediment studies; perform 
analyses and summarize results. 


G. R. Stephenson, Geologist Determine geologic and geomorphic 
parameters related to sediment yield. 


C. L. Hanson, Agr. Engin. Test various components in 
sediment models most applicable 
to rangelands. 


R. L. Engleman, Hydrol. Tech. Perform data. compilation and 
assist in analyses. 


Date of Initiation: September 1, 1969 


Expected Termination Date: Continuing 


INTRODUCTION 


Information on-sediment yield is almost entirely lacking for millions of 
acres of predominantly sagebrush rangeland under government land manage- 
ment and private ownership in the Northwestern United States. There is a 
growing concern for soil losses from intensively grazed rangelands, sedi- 
ment damage to reservoirs, and erosion of stream channels. 


Most rangeland watersheds in the Intermountain Northwest have large areas 
of relatively steep hillslope topography, and these areas need to be 
delineated for treatment to reduce erosion. Also, sediment yield infor- 
mation is needed for evaluating the benefits of watershed management and 
land treatment programs of the Bureau of Land Management and Soil Conser- 
vation Service. 


Range sites found in the Reynolds Creek Experimental Watershed represent 

a large percentage of the rangeland in the Northwest, and studies of sedi- 
ment yield are essential for developing sound management practices and 
planning appropriate multiple use of these lands. Good land management 
decisions require information on how vegetative changes, fencing, and 

land use alter the sediment yield potential of rangeland watersheds. The 
sources of these sediments need to be determined so that research data 

can be used to predict sediment yield for ungaged areas in terms of avail- 
able information on soils, climate, physiography, and use. Research is 
also needed to adapt the Universal Erosion Equation to rangelands. 
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Objectives: 


1. To determine the relationships between sediment yield and variables 
describing hydraulic and hydrologic factors and site and watershed 
characteristics that influence sediment yield. 


2. To formulate a sediment yield prediction procedure for rangeland 
watersheds in the Northwest. 


Suspended and bedload sediment yields from plots, channels, and water- 
sheds are measured by pumping sediment samplers, splitting devices, 
catchments, and hand sampling. Plots, microwatersheds, and watersheds are 
located on various soil types in different precipitation zones. A wide 
range of slope length, slope area, aspect, and relief ratio is represented. 
Rainfall intensity and duration data are available from a network of pre- 
cipitation gages, and snow data are available from snow courses, snow 
pillows, and other snow-measuring sites. Also, data on cover, 
topography, and soil factors, which influence erosion and sedimentation, 
are available through other investigations at the Northwest Watershed 
Research Center. 


PROGRESS 


Facilities and instrumentation for measuring runoff and sedimentation 

from six plots and two microwatersheds were completed in 1973. Descriptive 
data for the plot and microwatershed sites are contained in Table 1, and 
locations are shown on Figure 1 of the Introduction section of this report. 


if 
sary 


Precipitation at the Reynolds weather station during the 1973 water year 
was slightly above the 12-year average, as shown in Table 2. Monthly 
amounts were below or near normal except in April when the total was 

3.22 inches, the greatest monthly amount reported anywhere in Idaho. 

This record April precipitation also produced the peak runoff rates 

during the 1973 water year. Flooding did not occur since the precipi- 
tation was partly snow and there was no-frost in the soil. Consequently, 
sediment yields were low from most watersheds. Monthly runoff in 1973 

and the ll-year averages are compared in Table 2. Generally, near-average 
precipitation produced only about 60 percent of the average runoff. 


A study of winter storm and flood events on the Reynolds Creek Watershed 
shows that, during the 1962-72 period of record, most major floods occurred 
during December through March. Also, shorter-duration sedimentation 
records indicate that these events cause a large proportion of annual 
sediment yield on a water year basis. However, results of the study show 
the extreme variability in watershed areas that contribute runoff and sedi- 
ment during a single event. During the winter of 1973, secondary peak 
runoff rates occurred in mid-January, but these peak rates were lower 

than usual and runoff-erosion plots contributed no measurable sediment. 





OS Sa ea IRE AS AE OEY EA EET OTE AOE SC TL AEA ES GA SE er Ey 


Sk 
a eet 


ee en ne eee cae ne e  h 


44 


Bedload samplers have been procured to complement the suspended sediment 
sampling program. This equipment will provide total sediment yield data 
at a few streamflow sites, especially during major flood events. 


Processing and analysis of sediment data continued through 1973, and 
results will be submitted for publication in 1974. 


SIGNIFICANT FINDINGS 


The Reynolds Creek Experimental Watershed is well located to represent 
extensive rangeland areas in southwest Idaho and surrounding states to 
provide data on flood and sediment-producing events. Analyses of a few 
winter flood events have shown that runoff from portions of the 90-square- 
mile Reynolds Creek Watershed ranged from 0 to 0.8 inch. Similarly, 
measured suspended sediment concentrations ranged from 0 to 80,000 mg/1. 
The meteorologic, hydrologic, and hydraulic processes that cause runoff, 
erosion, and sediment movement are very complex and require extensive 
watershed data for analyses. 


Runoff and sediment yield were very low in 1973 when compared with 
average water year amounts since 1962; however, processing and analysis 
of sediment data are not complete. 


WORK PLAN FOR FY 1975 


The runoff and sediment collection program will continue at plots and 
watersheds, and additional sediment-measuring facilities will be con- 
structed at existing weirs. 


Bedload and suspended sediment samplers will be used to determine rates 
of sediment movement at selected sites during flood events. 


Sediment yield models will be investigated for applicability to Reynolds 
Creek Watershed conditions. 


———— 


—_eo 


REPORTS AND PUBLICATIONS 


Johnson, Clifton W., and Robin P. McArthur. 1973 
Winter storm and flood analysis, Northwest Interior. Proceedings 
Hydraul. Div. Spec. Conf., Am. Soc. of Civil Engin., Bozeman, MT. 
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TABLE 1.--Descriptions of runoff-erosion plots and microwatersheds, Reynolds 
Creek Watershed, 1973. (See Figure 1, Introduction, for locations.) 


No. . Elevation 


Summit 

Windy Point 
Nancy 
Nettleton 
Upper Sheep 
Reynolds Mtn. 


Microwatersheds 


Flats 
Nancy 





1/ 1972 Measurements 


TABLE 2.--Precipitation and runoff at selected stations, Reynolds Creek 
Watershed, 1973 water year. 


Precipitation~! Runoff = 


Outlet Weird! Reynolds Mtn.2/ 


1973 11-Year 1973 
Water Year | Average | Water Year 
(In) 

ral 

26 

; 233 

Ree | 

019 

20 

2.90 
75.33 

AG i 

14 

E(B be 

.05 

















12-Year 
Average 






Water 
Year 11.20 L135 3.08 eGo 20.36 13.24 


1/ Reynolds Weather Station, 076059, Weather Service Standard Gage 

2/ Reynolds Creek Outlet Weir, Station 036068, Watershed Area 90 sq. mi., 
Weir Elevation 3600 ft. 

3/ Reynolds Mountain Watershed, Station 166076, Drainage Area 100 a., F 
Weir Elevation 6600 ft. € 
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aN a) WATERSHED MODELING 


t 


Title: Developing, testing and evaluating watershed models. 
Personnel Involved: 
W. J. Rawls, Hydrol. Coordinate all modeling efforts 


and develop and field test water- 
shed models. 


C. L. Hanson, Agr. Engin. Develop and field test watershed 
models. 
G. R. Stephenson, Geologist Coordinate field data and test 


subsurface flow models 


sate Zuzel, Hydrol. Tech. Develop and field test watershed 
models. 


Date of Initiation: June 1973 


Expected Termination Date: Continuing 


e INTRODUCTION 


The rapid development of high-speed computers over the past decade has 
provided the impetus for significant advances in the development of models 
for synthesizing the hydrologic cycle. In the past few years, several 
notable contributions have appeared in the literature, such as Crawford 
and Linsley (1966); James (1970); Dawdy and O'Donnel (1965); Holtan and 
Lopez (1971); Sittner, Schauss, and Monro (1968); and TVA (1972). The 
models consist of mathematical expressions that define relationships and 
interactions among components of the hydrologic cycle. These models are 
operated with precipitation as input, and the model parameters are cali- 
brated to achieve an acceptable fit of simulated streamflow to observed 
streamflow. Different model parameter values are capable of yielding 
comparable simulations because many of the hydrologic processes are 
interrelated. Also, the various components may interact within the 
model, causing the parameters to compensate for each other and produce 
similar predictions. When such conditions exist, there is no positive 
verification that the mathematical expressions are a correct representa- 
tion of the various hydrologic processes. To develop better models, the 
important components need to be isolated and studied. 


Objectives: 


1. To make a comprehensive review of the existing watershed models. 
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2. To calibrate the most promising models, using Reynolds Creek V4 
Watershed data. - 


3. To test the sensitivity of the data inputs to the various models. 
4. To test the sensitivity of the various hydrologic components. 


5. To improve the representation of the hydrologic components found 
most important for Reynolds Creek data. : 


6. To use the models to study the influence of land use and con- 
servation practices on the hydrologic response of study watersheds. 


7. To apply the models to ungaged watersheds in the Northwest. 


PROGRESS 


The primary effort at the Northwest Watershed Research Center has been 
directed toward studying the individual components of the hydrologic 
cycle. Now, two types of hydrologic models are being studied to deter- 
mine the effect of each component. First, the distributed models, which 
chronologically route flows vertically and horizontally through the 
watershed; and second, the lumped models, which consider the watershed 
as a unit of hydrologic response. A comprehensive review of existing 
distributed watershed models applicable to the Northwest is in progress. és 
When this review is completed, the models will be evaluated according to «< 
their applicability to answer management problems. During the past 

year two lumped models and one distributed model have been developed 

and/or tested. -The following is a description of the models. 


Long-Term Snowmelt-Runoff Model 


A long-term snow-storage runoff-volume forecast model was developed 
using the equation: 


yrzart dix, + boxg + ..-- Dax, 


where, 
y = March-July runoff 
a = fitting coefficient 
X1°°Xy = peak snow course accumulation at each snow course 


by. .b, = coefficient for each snow course 
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In the forecast model each "b" coefficient would be realistic and 
hydrologically sound if it represented a percentage of the total 
drainage area. It can then be seen that the terms b,x] .... byxp 
in the forecast equation are analagous to a weighted average scheme 
and can be expressed as: 


n 
yrat 2 (by) (x) 


The hydrologic significance of b; can be maintained by the additional 
requirements of 


n 
ba.2,0,and.b ..bspsol 
i iz] 1 - 


provided that runoff and snow water contents are expressed in the same 
units, such as centimeters or inches. The basis for this reasoning is 
that the snow water content at a snow course is considered to be a 
discrete sample of snow water content associated with an area within a 
drainage basin, and not simply an index of runoff volume. These 
assumptions put the forecast equation on a sound hydrologic basis since 
snow water contents are now expressed as depth-area functions represent- 
ing the average snow water storage and can be related directly to runoff 
volume for any forecast period. 


An optimization technique that minimizes the error of y for a forecast 
period by searching for the optimum combination of "b" coefficients and 
the fitting coefficient "a" was used to calibrate the model to data. 


The forecast procedure was used to generate forecast equations for three 
drainage basins located in Southwest Idaho: Tollgate drainage (54 km) 
of the Reynolds Creek Watershed, the Middle Fork of the Boise River 
(2150 km“), and the entire Boise River above Boise, Idaho (7000 km?) . 
The average error for all forecast periods is 10 percent for the three 
watersheds. Comparisons with the accuracy of published forecasts for 
these forecast periods may be invalid due to the use of different 
parameters, but it does appear that the proposed method, using snow 
course data only, produces more accurate forecasts. 


Short-Term Snowmelt-Runoff Model 


A short-term snowmelt-runoff model was developed at the Northwest Water- 
shed Research Center (Cox and Zuzel, 1973). The model is basically a 
lag cross correlation model for which regression coefficients are cal- 
culated. The required inputs to the model are daily snowmelt from a 


snowmelt collector and daily streamflow. ; 


The daily forecasts produced by the model are based on the calculated 
regression coefficients. The generation of unrealistically high 
correlation coefficients, caused by the accumulation of input parameters, 
is not considered to seriously affect the accuracy of the model. 
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This paces has been tested on the Tollgate drainage of Reynolds Creek 
(54 km“), as well as larger drainage areas, including the. South Fork 
of the Boise River (1,022 km) and the Middle Fork of the Boise River 


(1,335 km“). 


Model output is mean daily streamflow 24 hours in advance. Errors in 
daily forecasts have been about 10 percent, and results indicate that 
forecasts up to 3 days in advance can be made with reasonable accuracy. 


Also, if the total water content remaining in the snowpack is used as 
input, extended volume forecasts can be made. This procedure has been 
tested and incorporated into the model. 


Subsurface Flow Model 


An example of a distributed model already tested is the transient, 
saturated-unsaturated, subsurface flow model developed by Freeze, 

1972. The model was tested in simulating the subsurface flow regime 

from snowmelt in a two-dimensional, vertical cross section on an 

instrumented slope in a small, upstream source area of Reynolds Creek. 

Field measurements show that the mechanism of streamflow generation is 
subsurface delivery of melt water over limited distances through shallow, 
high-permeability, low-porosity formations of altered and fractured € 
basalt. Overland flow is seldom observed on these slopes during ee 
snowmelt. Figure 1 shows the location of the study site, and Figure 2 

shows the hillslope topography, geology, and instrumentation. 


Output from the model gives plots of the pressure head, total head, 
and moisture content fields at any time step. The transient position 
of the water table and the development of the seepage face at the 
channel are determined from the pressure head field. Outflow through 
the seepage face is calculated from the head field. 


A simulation of the subsurface flow on this instrumented slope was 
run over a 100-day period during the 1971 spring snowmelt season. 
Figure 3 shows the comparison between the measured and simulated 
system for this run. The water levels in three piezometers and the 
estimated subsurface outflow to the stream, shown in Figure 3, are 
the comparative checks between field measurements and model output. 
The graphic comparisons are probably representative of the type of 
results that can be expected from the application of complex theoret- 
ical models to complex field situations under the constraints created 
by data limitations and the resulting uncertainties in the calibration 
procedures. 


eae yen 


a Ee a nN a ae a nh a mem a Oe an a eo ba nee 
whe —- hen ~ 


50 


References: 
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Forecasting runoff from universal surface gage snowmelt measure- 
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Crawford, N. H., and R. K. Linsley 1966 
Digital simulation in hydrology: Stanford Watershed Model IV, 
Tech. Report No. 39, Stanford Univ. Dept. of Civ. Engin. 


Dawdy, D. R., and T. O'Donnel 1965 
Mathematical models of catchment behavior, Proc. Jour. of the 
Hydraul. Div., ASCE, 91(HY4): 123-137. 


Freeze, R. Aw 1972 
Role of subsurface flow in generating surface runoff, 2. Up- 
stream source areas, Water Resources Res. 8(5): 1272-1283. 


Holtan, H. N., and N. C. Lopez 1971 
USDAHL-70 Model of watershed hydrology, Tech. Bul. No. 1435, 
ARS, USDA. 


James, L. Douglas 1970 
An evaluation of relationships between streamflow patterns and 
watershed characteristics through the use of Opset, a seif- 
calibrating version of the Stanford Watershed Model. Res. 
Report No. 36, Kentucky Water Resources Res. Inst., Lexington, 
Ree DOs 


Sittner, W. T., C. E. Schauss, and J. C. Monro 1968 
Continuous hydrograph simulation with an API type hydrologic 
model. Pres. at 49th Ann. Meeting, Amer. Geophys. Union, 
Washington, DC. © 


Tennessee Valley Authority 1972 
A continuous daily-streamflow model, Res. Paper No. 8, Knoxville, 
TN. 


SIGNIFICANT FINDINGS 


Long-term and short-term snowmelt forecast models were developed 
and successfully tested for a number of watersheds. 


A two-dimensional, transient model was used to simulate subsurface 
flow from a melting snowdrift to a stream channel. The results, when 
compared to the actual 1971 snowmelt season field measurements, were 
satisfactory. 
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WORK PLAN FOR FY 1975 

Snowmelt forecast models will be further refined. Also, the compre- ¢ 
hensive review of existing watershed models will be completed and the 
most promising will be calibrated to Reynolds Creek data. 


The subsurface flow model will be tested using the 1972 snowmelt data 
from Upper Sheep Creek. 


\ 


REPORTS AND PUBLICATIONS 


Stephenson, G. R., and Freeze, R. A. 1974 
Mathematical Simulation of Subsurface Flow Contribution to 
Snowmelt Runoff, Reynolds Creek Watershed, Idaho. Water Resources 
Res. In Press. 
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Figure 1. Location maps, Upper Sheep Creek study basin, Idaho. 
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Figure 2. Geologic cross section and guide to instrumentation 
along section AB in Figure l. 
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Figure 3. Calibrated transient flow system for the period April 5 . 
to July 23, 1971- 
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FRAIL LAND STUDIES 


Title: Hydrology of frail land watersheds 
Personnel Involved: 


C. W. Johnson, Supvy. Hydrol. Engin. Perform water balance and any 
other analysis needed. 


W. J. Rawls, Hydrologist Perform water balance and other 
analysis needed. 


Date of Initiation: October 1968 
Expected Termination Date: Field Work -— October 1973 


Interpretation and Summary - July 1974 


INTRODUCTION 


Proper management of watersheds for multiple use and protection requires a 
complete knowledge of the water balance with independent assessments of the 
water balance components. In order to determine the individual components 
for frail lands, a study of the hydrology of two watersheds, Rabbit Creek 
and Little Rabbit Creek, adjoining the Reynolds Creek Experimental Watershed, 
was initiated in 1968. A complete description of the watersheds and instru- 
mentation was included in the 1968-69 water-year report. 


Objective: 


To determine the water balance of two frail land watersheds in the North- 
west. 


PROGRESS 


Operation of five dual rain gage installations, two runoff stations, and 

10 soil moisture measuring sites were continued through October VoL iss to 
complete the 5-year study period. After October 1, 1973, all instrumentation 
was removed from the watersheds. The data for the 5-year period of study 
(1968-1973) has been error checked and processed. The termination report 

on the Rabbit Creek watersheds will be submitted to BLM about June 30, 1974, 
and will include: 


cS) 
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Summaries and Evaluation: <q 


as Evaluation and summary of the data collected. 
Zs Data and analysis of all storms causing runoff greater than 0.3 cfs. 
ous Frequency and intensity analysis of the precipitation data. 


4. A water balance for each watershed on a water-year basis. 


SIGNIFICANT FINDINGS 


None. 

WORK PLAN FOR FY 1975 
None. 

REPORTS AND PUBLICATIONS 
None. 
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